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Abstract 

The  Gurney  energy  method,  an  analytical  method  with 
which  to  compute  fragment  velocity  distributions  for  conti¬ 
nuous  casing  warheads,  is  reviewed  and  a  method  to  model 
velocity  losses  due  to  end  effects  on  cylindrical  charges 
and  end  projector  type  warheads  is  presented.  The  Taylor 
formula,  which  estimates  the  angle  of  projection  of  the 
fragments,  and  a  modified  Taylor  formula,  which  gives  better 
estimates,  are  also  presented.  The  method  is  extended  for 
use  with  casings  made  of  pre- formed  fragments.  Finally, 
four  conceptual  preliminary  designs,  based  on  information 
acquired  in  the  project,  are  investigated. 

The  modifications  made  to  the  equations  provide 
Improved  results  and  the  examples  confirm  that  the  modified 
Gurney  method  is  a  quick,  inexpensive  tool  for  use  in  pre¬ 
liminary  warhead  design. 


ESTIMATION  OF  VELOCITY  DISTRIBUTION 


OF  FRAGMENTING  WARHEADS  USING  A 
MODIFIED  GURNEY  METHOD 

I  Introduction 

Warheads  form  an  important  part  of  a  weapon  system  and 
proper  design  is  necessary  to  ensure  a  high  probability  of 
kill.  When  a  weapon  system  is  designed,  the  warhead  is 
generally  not  an  "off-the-shelf"  item  and  it  is  designed 
with  the  intent  of  defeating  a  particular  class  of  targets. 
The  damage  mechanisms  of  conventional  warheads  are  blast  and 
fragments.  In  some  cases,  such  as  air-to-air  missiles,  it 
is  the  fragments  generated  from  the  explosion  of  the  war¬ 
head  which  will  incapacitate  or  destroy  the  target  by 
piercing  or  passing  through  it.  The  destructive  power  of 
the  fragments  is  actually  derived  from  their  kinetic  energy 
and,  in  order  to  estimate  that  energy,  the  mass  and  velocity 
of  the  fragments  must  be  known.  The  velocity  and  mass 
distribution  of  fragments  from  exploding  warheads  is  a  com¬ 
plex  function.  This  is  exemplified  in  Fig  1  where  a  20-mm 
round  is  loaded  with  tetryl  and  is  detonated.  The  flash 
radiographies ,  which  have  been  touched  up  to  ease  reproduc¬ 
tion,  show  that  the  round  expands  circumferentially  and 


1 


Fig  1  (Cont'd).  Detonation  of  Tetryl  Filled  20  mm  Projectile  (Ref  18  ) 
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eventually  breaks  into  many  fragments  of  differing  shapes, 
and  masses.  This  is  the  typical  behavior  of  an  exploding 
warhead. 

Of  prime  importance  to  the  designer  is  the  velocity 
distribution,  which  he  must  be  able  to  estimate  to  a  reason¬ 
able  degree  of  accuracy  to  determine  the  damage  that  the 
warhead  can  inflict.  There  are  essentially  three  methods 
which  can  be  used  for  warhead  design.  The  first  is  the 
experimental  method  where  a  warhead  of  a  certain  design  is 
tested  and  modified  until  the  desired  characteristics  are 
attained.  A  second  method  is  to  use  computer-aided  design. 
Several  computer  programs  exist  which  estimate  velocity  dis¬ 
tributions  quite  well;  an  example  is  the  HEMP  code  (Ref  9). 
Unfortunately,  both  methods  are  time  consuming  and  can  be 
very  expensive.  The  third  method  is  a  theoretical  method 
which  is  relatively  simple,  accurate  for  cases  where  L/D>2, 
and  is  well  suited  to  preliminary  design  of  warheads.  It  is 
known  as  the  Gurney  energy  method.  Because  of  the  simplis¬ 
tic  approach  taken  in  deriving  the  relationship  between 
fragment  velocities,  Gurney  energy  and  C/M,  the  Gurney 
energy  method  holds  but  for  . l^C/M^S  (Ref  6,13).  For  values 
of  C/M  outside  these  limits,  gas  dynamic  equations  must  be 
used  to  compute  reasonable  velocity  estimates.  However,  most 
warheads  used  in  modern  weapons  fall  within  that  range  of 
C/M's  and  the  Gurney  method  can  be  used  to  predict  fragment 
velocity  distributions,  which,  in  most  cases,  are  within  ten 
percent  of  experimental  data.  However,  experimental  data  are 
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^  difficult  to  obtain  accurately  because  of  the  nature  of  the 

process.  The  fragments  move  so  rapidly  that  it  is  nearly 
impossible  to  measure  exact  velocities.  Also  in  cases  where 
length  to  diameter  (L/D)  of  the  warhead  is  less  or  equal  to 
two,  estimates  of  the  velocity  distributions  calculated  by 
Gurney's  method  are  in  error  by  as  much  as  30  percent  at 
the  free  ends. 


r> 


The  purpose  of  this  study  is  to  collect  and  present 
currently  available  information  on  the  Gurney  method  for 
estimating  velocity  distributions  for  warheads  with  contin¬ 
uous  casings  and  to  extend  it  for  use  on  warheads  made  of 
pre-formed,  or  discrete,  fragments.  Also,  the  equations 
will  be  modified  to  account  for  variations  in  velocities  due 
to  end  effects,  which  are  especially  important  in  warheads 
where  L/DS2.  Finally,  the  material  developed  and  presented 
will  be  used  in  the  conceptual  preliminary  design  of  three 
example  warheads  to  illustrate  its  application  to  the  design 
of  different  warhead  geometries.  This  improved  method 
should  be  used  for  preliminary  design  because  it  is  quicker, 
easier  to  use  and  cheaper  than  the  other  methods. 

In  addition  to  the  Gurney  method,  Taylor's  formula, 
which  estimates  the  angle  of  departure  of  the  fragments,  will 
be  reviewed.  A  more  accurate  equation  which  estimates  the 
angle  of  departure  will  be  presented. 
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Modification  to  Gurney  and  Taylor  Equations 

The  modifications  to  the  Gurney  method  are  derived  from 
available  experimental  data.  From  that  data,  the  Influence 
of  end  effects  on  the  velocity  distribution  of  the  fragments 
from  cylindrical  and  end  projector  type  warheads  can  be 
observed.  Terms  in  the  Gurney  equation  for  a  particular 
simple  metal/explosive  geometry  will  be  modified  such  that 
the  computed  velocity  distribution  approximates  the  test 
data  for  that  particular  geometry.  These  data  are  contained 
in  Ref  2,  9,  10,  11  and  19.  Although  there  is  not  an  abun¬ 
dance  of  experimental  data,  it  is  sufficient  to  determine 
appropriate  modifications  to  the  Gurney  equations .  The 
modified  equations  are  Intended  for  use  in  preliminary 
design,  after  which  experimental  testing  or  computer  simu¬ 
lation  would  be  used  to  finalize  the  design.  Some  conceptual 
preliminary  designs  of  end  projector  type  warheads  will  be 
investigated.  This  type  of  warhead  is  particularly  suited 
to  air-to-air  missiles  using  proportional  navigation.  Being 
a  directional  warhead,  its  fragments  would  travel  approxi¬ 
mately  along  the  line  of  sight  to  the  target  and,  if  detonated 
close  enough  to  the  target,  should  give  higher  probabilities 
of  kill  than  a  similar  missile  using  an  isotropic  warhead. 
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II  Review  of  the  Derivation  of  Gurney  and  Taylor  Formulae 


Gurney  Equations  (Ref  6,  Ref  13) 

The  following  developments  may  be  found  in  the  cited 
references  and  are  presented  here  for  the  sake  of  complete¬ 
ness  and  convenience  of  the  reader. 

The  Gurney  method  is  straightforward  and  is  based  on 
energy  and  momentum  balances.  It  is  accurate  over  a  wide 
range  of  charge  to  mass  ratios  (C/M)  and  works  best  for  one¬ 
dimensional  translation  of  a  metal  surface  (Ref  10).  The 
method  makes  a  number  of  assumptions  which  simplify  the 
problem.  First,  in  a  metal/explosive  system,  the  chemical 
energy  of  the  explosive  is  assumed  to  be  completely  con¬ 
verted  to  kinetic  energy  upon  detonation  of  the  explosive. 
This  energy  imparts  a  velocity  to  the  product  gases  and  the 
metal  casing.  Secondly,  the  velocity  profile  of  the  product 
gases  is  linear  and  it  is  constant  throughout  the  metal 
thickness.  Representation  of  this  is  shown  in  Fig  2. 
Thirdly,  charge  gases  after  detonation  are  assumed  to  be 
equally  dense  everywhere  and  are  expanding  uniformly. 

Lastly,  rarefaction  waves  which  are  created  behind  the  reac¬ 
tion  zone  are  neglected. 

The  assumptions  made  enable  us  to  reduce  what  would  be 
complex  equations  to  simple  ones  which  can  be  solved  to 
yield  reasonable  results  for  velocity  distributions. 

Although  the  chemical  energy  Is  actually  transformed  to 
kinetic  and  thermal  energy  and  light,  the  kinetic  energy 
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Fig  2.  Linear  Gas  Velocity  Profile  in  Metal/Exploaive  System 


term  greatly  overpowers  all  other  terms,  except  for  nuclear 
explosions,  which  will  not  be  considered.  The  assumption 
about  constant  density  is  far  from  reality  because  the  gases 
near  the  reaction  zone  will  be  denser.  Henry  (Ref  6) 
compared  the  Gurney  method  with  one  based  on  a  parabolic 
distribution  of  charge  gas  density.  Results  showed  that  the 
curves  for  characteristic  velocity  of  the  fragments  coincide 
for  values  of  C/M  up  to  1.0,  whereas  at  C/M  of  5.0,  the 
difference  is  approximately  12  percent.  However,  experi¬ 
mental  data  are  not  available  to  support  or  contradict 
Henry's  results.  The  added  complexity  of  assuming  a  para¬ 
bolic  distribution  is  not  worthwhile  and  a  constant  density 
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assumption  permits  an  easy  derivation  of  a  simple  relation¬ 
ship  between  velocity,  Gurney  energy  and  C/M,  which  gives  a 
reasonably  accurate  fragment  velocity  distribution. 

Finally,  the  last  assumption  about  rarefaction  waves  is  also 
misleading.  Upon  initiation  of  the  explosive  by  some  ade¬ 
quate  means,  a  detonation  wave  travels  into  the  explosive, 
and  rarefaction  waves  travel  in  the  opposite  direction 
(Ref  22).  The  effect  of  the  rarefaction  waves  is  to  de¬ 
crease  the  pressure  which  builds  up  to  accelerate  and 
rupture  the  casing.  Although  the  last  two  assumptions  do 
not  truly  represent  reality,  they  constitute  an  important 
part  of  Gurney's  method  by  introducing  two  cancelling 
errors;  first,  assuming  constant  density  causes  the  velocity 
estimate  to  be  low;  secondly,  neglecting  the  existence  of 
rarefaction  waves  causes  the  velocity  estimate  to  be  high. 
This  results  in  the  method  being  accurate  over  a  wide  range 
of  C/M  ratios  (0.1  to  5.0) (Ref  6,13). 

The  Gurney  method  determines  the  initial  velocity  of 
the  metal  casing  as  a  function  of  C/M  and  the  Gurney  spe¬ 
cific  energy,  JTE.  The  functions  differ  for  various  warhead 
geometries,  but  only  involve  C/M  and  S2E.  The  quantity 
/ZE  is  determined  empirically  and  is  a  characteristic  of 
each  explosive.  It  has  been,  measured  for  several  commonly 
used  explosives.  The  values  for  the  Gurney  energy  quoted  in 
this  report  were  obtained  from  Ref  A  and  13.  For  some 
explosives,  two  or  more  different  values  are  given,  illus¬ 
trating  the  difficulty  encountered  in  obtaining  exact  data. 
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However,  the  difference  between  the  values  Is  slight.  When¬ 
ever  Gurney  energy  values  were  required  for  computations, 
the  smallest  values  found  In  the  references  were  used 
because  the  Gurney  method  overestimates  the  velocity  in  most 
Instances  and  by  using  the  smallest  value  of  /ZE,  agreement 
between  experimental  data  and  the  Gurney  estimate  was  better. 
In  analyzing  the  experimental  data  used  later  in  this  report, 
it  was  necessary  to  estimate  the  value  of  the  Gurney  energy 
of  explosives  for  which  no  data  were  found.  Since  the 
energy  appears  to  vary  linearly  with  the  detonation  velocity 
of  the  explosive,  a  linear  least  squares  fit  was  used  to 
obtain  the  energy  values.  The  graph  is  shown  in  Appendix  A. 
The  experimental  data  were  available  for  cylindrical  charges 
and  sandwich- type  configurations;  consequently,  derivation 
of  the  Gurney  equations  for  those  geometries  will  be 
reviewed.  Equations  for  other  geometries  can  be  found  in 
Ref  6  and  13. 

A  cylindrical  charge  is  illustrated  in  Fig  3.  The 
cylinder  is  assumed  infinite  in  length  so  that  end  effects 
can  be  neglected. 

From  the  conservation  of  energy  principle,  the  total 
chemical  energy  of  the  explosive  before  detonation,  CE,  is 
approximately  equal  to  the  kinetic  energy  after  detonation 
(thermal  and  light  energy  are  neglected) . 

Chemical  Energy  m  Kinetic  Energy  ,  Kinetic  Energy 

of  Explosive  of  Metal  of  Product  Gases  (1) 


Fig  3.  Cylindrical  Casing 

R 

CE  =  %MV2  +  %Jv*(r)27rrp(r)dr  (2 

o 

C  “  mass  of  explosive  per  unit  area 

E  *  energy  per  unit  mass  of  explosive 

M  *  mass  of  metal  per  unit  area 

V  ■  initial  velocity  of  metal  imparted  by  the 
explosion 

)  “  density  of  the  charge  gas  at  some  point  r 

R  -  displacement  of  the  casing  at  some  time  after 
detonation 

r  *  distance  from  centerline 


Hence,  for  a  cylinder  of  unit  length,  the  mass  of  the 
charge  gas  between  r  and  r+dr  Is  2irrp(r)dr.  The  density, 
assumed  constant,  at  r  is  given  by 


p(r)  -  C/nR* 


(3) 


and  the  velocity  at  r  is 


v(r) 


V 


(4) 


O 


because  of  the  assumptions  of  constant  density  and  linear 
velocity  profile. 

Substituting  for  p(r)  and  v(r)  into  Eq  (2)  and  inte¬ 
grating,  we  obtain 


CV* 

CE  -  %MVJ  + 


(5) 


from  which 


M 


S2E 


H  +  I 


(6) 


or 


( 


/2E  [ _ C/M _ ]i 

[l  +  .5  C/MJ 


(7) 
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N 


Fig  4,  Asymmetric  Sandwich 

C 

The  equation  for  an  asymmetric  sandwich  is  developed 
as  follows: 

Figure  4  represents  an  asymmetric  sandwich.  The  sand- 
^  *wich  is  assumed  infinite  in  length  so  that  end  effects  can 
be  neglected. 

After  detonation,  there  will  exist  a  plane  where  the 
gas  is  stationary;  let  this  plane  be  at  a  distance  a  from 
the  surface  of  the  plate  of  mass  M  per  unit  area  and  at  a 
distance  b  from  the  surface  of  the  plate  of  mass  N  per  unit 
area.  Applying  the  energy  conservation  principle  to  before 
and  after  the  detonation, 

C 
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Chemical  Energy 
of  Explosive 


Kinetic  Energies 
of  Plates  M  and  N 


Kinetic  Energies  of 
+  Charge  Gases  Between 
Plates  M  and  N 


(8) 


that  is  , 

a  b 

CE  -  %MV*  +  %Jvl (y)p(y)dy  +  %NV’  +  %Jv* (y)p(y)dy  (9) 

o  o 


where 

VM  «  initial  velocity  of  plate  of  mass  M 
VN  «  Initial  velocity  of  plate  of  mass  N 

The  remaining  terms  are  defined  as  in  Eq  2  and  primes  indi¬ 
cate  the  negative  y  direction. 

Now  from  the  conservation  of  momentum. 


Momentum  of  +  Momentum  of  gases 
Plate  M  between  o  and  a 


Momentum  of 
Plate  N 


Momentum  of 
+  gases  between 
o  and  b 


(10) 


that  is , 


a  b 

MVm  +  |p(y)v(y)dy  -  NV^  +  jp(y)v(y)dy 
o  o 


(ID 


but  since  the  gases  are  assumed  to  expand  uniformly 
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r 


v  ■  — 

VM  b 


which  implies  that 


v(y)  -  (I]VM 

**>  ■  (f]v» 


(12) 


(13) 


(14) 


Also,  the  density  was  assumed  to  be  constant,  implying  that, 


p(y)  ■  p(y)  ■ 


(a+b) 


(15) 


Substituting  Eqs  (12),  (13),  (14)  and  (15)  into  (11)  and 
integrating  yields 


MV  +  — -  V  -  NV 
M  a+b  2  "  M 


C  b2 
a+b  ~T~  a 


(16) 


solving  for  b/a  gives 


b  _  C/M  +  2 
a  C/M  +  2N/M 


(17) 


Substituting  Eqs  (13),  (14)  and  (15)  into  Eq  (9)  and  solving 
(  gives 
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2CE  „  M  +  JL.  £  +  N  ^.  +  JL.  ^  (18) 


w 


a+b  3 


a+b  3aJ 


Solving  for  VM  gives 


-  vTE 

rl 


1  l-Kb/a)* 
3  l+(b/a) 


Nfk]2  +  M 
cUi  c 


-*5 


(19) 


n 


The  value  of  VN  is  determined  from  Eq  (12) . 

If  Eq  (17)  is  substituted  into  Eq  (19),  we  obtain 


where 


v„  -  M 


1  +  A1  +  »  Az  +  M 


3  ( 1+A)  C 


(20) 


C/M  4*  2  _  b 


C/M  +  2N/M  a 


(21) 


The  equation  for  an  open-faced  sandwich  can  be  easily 
obtained  by  letting  N-0  .  This  gives 


/Z  (C/M) 


V  -  /ZE  - - 

M  ((C/M)2  +  5  (C/M)  +  4)*s 


(22) 


We  have  thus  derived  the  Gurney  equations  which  will  be 
necessary  to  analyze  the  experimental  data. 
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As  was  stated  earlier,  Taylor's  formula  predicts  the 
angle  of  departure  of  the  fragments.  This  formula  is 
derived  as  follows .  Figure  5  shows  the  detonation  of  a 
charge  against  a  metal  plate.  The  detonation  velocity,  D, 
is  parallel  to  the  surface  of  the  plate.  It  is  assumed 
that  acceleration  of  the  plate  to  its  final  velocity  is 
instantaneous  and  steady  state  is  achieved.  The  plate  is 
assumed  to  undergo  no  net  shear  flow  and  therefore  there  is 
no  change  in  length  or  thickness.  Although  the  plate 
accelerates  gradually  to  its  final  velocity,  VM,  the  time 
lapse  is  so  short  that  for  practical  purposes  it  can  be 


Fig  5.  Projection  of  Metal  Plate  by  Detonation  of  Explosive 


assumed  to  be  instantaneous.  At  some  time  after  detonation, 
the  deflection  of  the  plate  from  its  original  position  is 
a  degrees.  The  point  B  therefore  goes  to  point  B  and 


AB 


(23) 


The  line  AC  is  drawn  perpendicular  to  BB.  Therefore, 

*  * 

angle  CAB  »  angle  CAB  -  a/2,  since  ABB  is  an  isosceles 

triangle.  The  velocity  vector  D  is  along  AB  and  the  velocity 

> 

vector  Vj.  is  along  BB.  From  trigonometric  and  geometric 
relations,  we  obtain 


V  /2 

Sin  a/2  -  -5—  (24) 

O  ” 

which  is  Taylor's  formula. 

The  velocity  VM  can  be  obtained  from  Gurney's  equation 
for  the  appropriate  geometry. 

Accuracy  of  the  Equations 

As  stated  before,  the  Gurney  equations  are  reasonably 
accurate  for  a  large  number  of  configurations ,  L/D  and  C/M. 
In  a  cylindrical  bomb  of  constant  diameter  and  L/D>2,  the 
initial  fragment  velocity  is  essentially  the  same  for  all 
the  fragments  in  the  central  part  of  the  casing,  neglecting 
the  small  portion  of  fragments  generated  at  the  ends.  How¬ 
ever,  in  many  cases,  the  diameter  of  the  explosive  charge 


§. 
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is  not  constant.  Many  projectiles  have  ogival  noses.  It 
is  therefore  more  appropriate  to  take  a  local  charge  to 
mass  ratio  for  use  in  Gurney's  equations.  These  equations 
were  derived  using  an  element  of  explooive/metal  system  of 
unspecified  size.  For  example,  to  use  the  total  C/M  for  a 
cone  shaped  warhead  to  find  the  velocity  distribution  of 
the  whole  casing  would  be  erroneous.  The  fragments  at  the 
apex  of  the  cone  would  have  a  much  lesser  velocity  than  the 
ones  at  the  base,  but  using  the  total  C/M  would  predict  the 
same  velocities  for  all  the  fragments.  Consequently,  care 
must  be  taken  in  applying  Gurney’s  formula  to  ensure  that 
the  correct  C/M  is  used.  The  accuracy  of  the  method  can  be 
observed  in  Fig  6  where  the  velocity  distribution  obtained 
by  Gurney's  method,  using  local  C/M  along  the  projectile,  is 
compared  to  experimental  data  for  a  typical  high  explosive 
projectile  (Ref  10,11).  As  can  be  seen,  the  agreement  is 
quite  good.  It  must  be  noted  that  in  this  case  the  length 
to  diameter  ratio  is  fairly  large. 

At  either  end,  where  the  C/M  is  small,  Gurney's  equa¬ 
tion  predicts  a  lower  value  than  at  the  center  of  the 
projectile  and  this  is  borne  out  by  the  data.  It  is  expected 
that  similar  agreement  should  occur  for  other  warheads  of 
similar  configurations.  Unfortunately,  it  has  been  observed 
that  Gurney's  method  does  not  work  well  for  configurations 
where  L/D-2,  even  if  local  C/M's  are  used.  The  estimates  of 
the  fragment  velocities  at  the  ends  are  in  error  by  as  much 
as  28  percent  in  some  cases.  Figure  7  compares  the  estimate 
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from  Gurney* 8  equation,  Eq  (7),  with  experimental  data  for 
a  cylindrical  charge  of  L/D-2.  Although  the  C/M  is  the  same 
for  the  whole  length  of  the  cylinder,  which  implies  that  the 
velocity  should  be  the  same  for  all  the  fragments ,  it  can  be 
seen  that  the  fragment  velocities  at  the  ends  are  less  than 
at  the  center.  The  Gurney  equation  provides  a  reasonable 
estimate  for  fragments  located  from  about  0.5  to  0.8  of  the 
relative  distance  from  the  initiated  end,  but  fails  for  the 
rest  of  the  casing.  Taylor's  formula  also  shows  agreement 
and  large  errors  within  roughly  the  same  areas,  which  is  to 
be  expected,  since  its  estimate  depends  on  the  velocity 
obtained  from  Gurney's  equation. 

The  equation  loses  its  accuracy  at  the  ends  because  of 
the  free  end  effects.  These  end  effects  result  from  rarefac¬ 
tion  waves,  which  are  neglected  in  Gurney's  theory,  but  it 
appears  that  they  must  be  taken  into  account.  As  the  Gurney 
equations  are  presented,  they  are  not  adequate  for  evalua¬ 
ting  velocity  distributions  for  warheads  with  L/DS2.  Com¬ 
puter  codes  are  more  accurate.  The  best  HEMP  code,  which 
assumes  elastic-plastic  properties  for  the  casing  and 
accounts  for  gas  leakage  between  the  fragments,  models  the 
experimental  data  quite  accurately  (Ref  9).  Unfortunately, 
it  is  a  complex  tool  requiring  much  computer  time  to  obtain 
the  results.  A  method  of  modifying  the  Gurney  equation  for 
cylindrical  charges  so  that  it  computes  more  accurate  results 
is  investigated  in  the  next  section.  Use  of  an  appropriately 
modified  Gurney  equation  would  save  time  and  money.  Pre- 
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liminary  design  could  be  done  with  this  equation  and  the 
final  design  confirmed  with  one  or  two  runs  of  the  HEMP 
code . 
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End  Effects 

End  effects  are  present  in  all  cases  but  have  a  varying 
degree  of  importance  in  computing  the  velocity  distribution 
of  the  fragments.  In  cylindrical  warheads  of  L/D>2,  the 
end  effects  could  be  neglected  and  the  error  in  computing 
the  velocity  distribution  would  be  small.  Most  of  the  frag¬ 
ments  in  such  cases  would  attain  the  maximum  velocity  which 
can  be  achieved  for  a  given  C/M  and  L/D.  The  fragments  at 
the  ends  would  have  reduced  velocities  but  since  L/D  is 
large,  the  overall  percentage  of  those  fragments  would  be 
small.  Gurney's  equation  would  give  a  reasonable  estimate 
of  the  velocity  distribution,  being  in  error  only  at  the 
ends.  However,  for  design  purposes,  it  is  desirable  to  be 
able  to  estimate  the  velocity  distribution  for  all  the  frag¬ 
ments.  In  cases  where  L/D^2,  the  influence  of  end  effects 
are  such  that  few  of  the  fragments  reach  velocities  pre¬ 
dicted  by  Gurney's  equation  for  a  cylindrical  charge.  Conse¬ 
quently.  a  method  which  yields  better  estimates  is  required. 

End  effects  result  from  the  rarefaction  waves  generated 
in  the  system.  When  detonation  occurs,  a  shock  wave  passes 
through  the  explosive,  transforming  the  chemical  energy  of 
the  explosive  to  kinetic  energy  and  a  reaction  zone  of 
finite  thickness  forms  behind  the  shock  wave.  The  pressure 
and  temperature  rise  sharply  behind  the  shock,  but  since 
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there  is  no  interface  to  maintain  this  rise,  they  begin  to 
fall  off.  The  explosive  process  continues  through  the 
reaction  zone  and  beyond  into  the  burnt  gases  so  that  rare¬ 
faction  waves  are  created  (Ref  21).  The  process  is  sketched 
in  Fig  8.  Because  of  the  high  detonation  velocities  in¬ 
volved  in  the  process,  the  explosive  matter  is  completely 
detonated  before  the  casing  breaks  apart.  The  rarefaction 
waves  formed  at  the  initiated  end  of  the  explosive  travel 
in  a  direction  opposite  that  of  the  detonation  shock.  The 
center  part  of  the  detonation  products  is  still  confined 
and  the  pressure  loss  is  less  there  than  at  the  initiated 
end.  When  the  detonation  shock  reaches  the  opposite  end, 
the  gases  are  no  longer  confined  and  the  rarefaction  waves 
cause  the  loss  of  pressure.  In  cases  where  L/D  is  small,  the 
rarefaction  waves  reach  farther  inside  the  casing  towards 
the  middle  and  the  pressure  is  decreased  over  a  wider  axial 
distance.  Because  the  pressure  is  decreased,  the  velocity 
of  the  fragments  will  be  less.  This  loss  in  velocity  is 
similar  to  having  a  smaller  C/M  at  the  ends,  that  is,  the 
velocity  is  proportional  to  the  C/M.  Therefore,  end  effects 
can  be  modelled  by  decreasing  the  local  C/H  along  the  casing 
by  an  appropriate  amount.  The  difficulty  lies  in  finding 
the  amount  of  the  decrease.  Whatever  the  modification 
required  is,  it  must  give  reasonable  results  and  be  of  a 
simple  form  so  that  it  is  easy  to  apply.  Returning  to  Fig  8, 
It  can  be  seen  that  the  pressure  loss,  consequently  the  loss 
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Shock  Wave 


a)  Progression  of  Detonation 


b)  Complete  Detonation  of  Explosive. 
Casing  has  not  burst  yet. 


Fig  8.  Explosive  Process 
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in  velocity,  varies  with  the  position  along  the  casing  and 

» 

is  greater  at  the  ends.  End- Initiated  cylinders  are  inves¬ 
tigated,  since  experimental  data  are  available  for  these. 

Modification  to  the  Gurney  Formula  (Cylinders) 

Reduction  of  the  local  C/M  can  be  achieved  by  assuming 
that,  for  confutations!  purposes,  an  amount  of  explosive  is 
removed  from  the  ends.  Because  simplicity  is  a  factor,  simple 
geometric  shapes,  such  as  cones  or  hemispheres,  were  investi¬ 
gated  as  likely  amounts  to  be  removed.  The  method  was  to 
assume  the  removal  of  a  cone  or  hemisphere  of  certain  dimen¬ 
sions  from  either  end,  confute  the  velocity  distribution 
using  Gurney's  formula  for  a  cylindrical  case,  with  the  appro¬ 
priate  local  C/M  and  compare  the  results  obtained  to  the  expe¬ 
rimental  data  for  that  particular  combination  of  explosive 
type,  L/D,  and  C/M.  Best  agreement  was  obtained  when  two 
cones,  one  of  height  2R  and  of  base  equal  to  the  diameter, 
and  one  of  height  R  and  of  base  equal  to  the  diameter,  were 
removed  from  the  initiating  end  and  the  opposite  end,  res¬ 
pectively  (Ref  20).  Figure  9  illustrates  this.  The  cylinder 
is  divided  into  three  regions.  In  region  B,  no  modification 
is  required  and  the  correction  factor  is  1.0.  In  region  A 
and  C,  the  factor  is  determined  as  follows. 

In  region  A, 


r  _  2R-x  _  ,  x 

TR  1  ~  7F  (25) 
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End  View 


Modification  for  Cylinder  Filled  with  Explosive 


If  the  cylinder  is  viewed  from  Che  end,  the  proportional 
amount  of  explosive  at  that  location  is 


F(x)  -  area  of  sector  S  _  ttR2  -  nr2 
total  area  nR* 


(26) 


Similarly  in  region  C,  the  amount  of  explosive  at  any 
one  point  would  be 


F(x)  -  1  -  1  - 


(27) 


Therefore,  the  correction  factor  for  the  whole  cylinder  can 
be  written  in  a  compact  form,  suitable  for  insertion  in  a 
computer  program,  as 


F(x)  -  1 


(28) 


Multiplying  C/M  by  F(x)  will  give  the  appropriate  local  C/M 
to  use  in  Gurney's  equation.  The  cone  at  the  initiating  end 
is  greater  because  it  allows-  for  the  build-up  of  the  explo¬ 
sion  at  that  end. 

The  Gurney  equation  to  be  used  for  cylindrical  cases 
when  L/D^2,  filled  with  any  explosive,  and  with  . 1SC/M55 
becomes 
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VM(x) 


m\ _ F<*)c/H _ I* 

ll  +  ,5F(x)C/M) 


(29) 


Modification  to  the  Taylor  Formula  (Cylinders)  (Ref  20) 

In  Taylor's  derivation  for  deflection  angle,  steady 
state  conditions  were  assumed  to  exist.  In  many  cases,  this 
is  incorrect  because  there  exists  a  velocity  gradient  in  the 
metal.  It  is  therefore  necessary  to  use  a  different  formula 
to  obtain  the  projection  angles.  The  equation  proposed  by 
Banders -Pherson  gives  reasonable  agreement  between  calcu¬ 
lated  values  and  experimental  data.  In  his  derivation,  he 
considered  two  closely  spaced  points  along  the  metal  surface. 
Each  point  was  assumed  to  accelerate  according  to  the 
equation 

v  -  <3o> 

where  V  is  the  velocity  at  any  instant,  Vp  is  the  ultimate 
velocity  which  the  element  would  attain,  T  is  the  time  at 
which  the  detonation  front  reaches  the  element,  and  t  is  the 
time  constant  of  acceleration,  i.e.  the  time  required  for 
the  element  to  reach  VQ (1-1/e).  V  ,  T,  and  t  all  vary  with 
initial  location.  If  there  is  no  velocity  gradient,  due  to 
the  acceleration  of  tha  metal,  Vo  and  t  are  independent  of 
initial  location  and  the  two  elements  follow  identical 
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trajectories.  That  equation  is 


V  V't 
Sin  a  -  J  -  1*1 


<V>‘ 


where 


•  angle  from  vertical  (see  Fig  6) 

•  fragment  velocity,  obtained  from  the  modified 
Gurney  equation 

-  detonation  velocity  of  the  explosive 

•  velocity  gradient  in  the  metal 

•  time  constant  of  acceleration 


The  equation  was  derived  empirically.  The  velocity  gradient 
V'  ,  is  the  difference  in  velocity  between  two  adjacent 
points,  divided  by  the  distance  separating  them,  and  the 
value  of  t,  defined  as  the  time  for  the  element  to  reach 
VQ(l-l/e),  is  estimated  from  computing  the  velocity  of  the 
fragment  at  different  times  using  experimental  data  for 
radial  expansion  of  cylindrical  charges  in  the  standard 
cylinder  test.  It  may  be  noticed  that  if  there  is  no  veloc¬ 
ity  gradient,  the  equation  becomes  Taylor’s  formula. 

Modification  for  Pre-Formed  Fragments  (Cylinders) 


From  explosive  tests  with  steel  cylinders  with  contin¬ 
uous  casings,  the  first  appearance  of  gas  leakage  is  at 
expansion  ratios  of  1.6  to  2.1,  depending  on  the  type  of 
explosive  filling,  i.e.  its  power  (Ref  10).  In  similar 
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tests  with  cylinders  made  of  pre-formed  fragments,  gas 
leakage  occurs  at  expansion  ratios  of  1.18  to  1.26  (Ref  19). 
This  Implies  that  much  of  the  energy  escapes  between  the  frag- 
ments  and  Is  rot  available  to  accelerate  the  fragments .  This 
can  be  represented  by  reducing  the  overall  C/M  for  the 
cylinder.  Therefore,  the  correction  was  obtained  by  reducing 
the  total  C/M  until  acceptable  agreement  between  computations 
and  experimental  data  was  obtained. 

Modification  to  the  Equation  for  End  Projectors 

One  of  the  purposes  of  this  study  is  to  apply  a  modi¬ 
fied  Gurney  method  to  the  design  of  end  projectors.  End 
projectors  of  Interest  are  of  circular  4hape  (to  fit  Inside 
an  air-to-air  missile)  and  are  made  of  pre-formed  fragments. 
Such  warheads  are  also  subject  to  end  effects;  however,  the 
modification  will  be  different  since  the  configuration  is 
different  from  that  of  the  cylindrical  case.  In  an  end  pro¬ 
jector,  the  detonation  wave  is  parallel  to  the  layer  of 
fragments  to  be  projected,  instead  of  being  perpendicular  as 
it  is  in  a  cylinder.  Figure  10  Illustrates  a  cross  section 
of  a  sample  end  projector.  Since  end  effects  were  modelled  by 
reducing  the  local  C/M  in  the  case  of  cylinders,  it  is  appro¬ 
priate  to  do  the  same  for  this  geometry,  because  the  presence 
of  rarefaction  waves  will  reduce  the  fragment  velocities 
at  the  edges  of  the  projector  like  they  did  at  the  ends  of 
cylinders.  In  this  case,  the  Gurney  equation  to  be  used  is 
that  for  an  open-faced  sandwich  (see  Fig  10).  Since  the 
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geometry  involved  is  different,  the  correction  for  end 
effects  must  be  slightly  different.  Again  the  reduction  of 
local  C/M  is  to  be  modelled  using  a  simple  geometric  shape 
It  was  found  that  best  agreement  with  experimental  data 
(Ref  11)  occurred  when  a  cone  of  explosive  of  height  D/2  and 
base  D  was  assumed  to  be  the  amount  of  explosive  present 
(See  Fig  10).  Also,  since  the  end  projector  is  made  of 
pre~ formed  fragments,  the  total  C/M  is  reduced  to  account 
for  gas  leakage  between  the  fragments.  The  reduction  is 
assumed  to  be  the  same  amount  as  that  for  cylindrical  case. 

To  recapitulate,  the  Influence  of  end  effects  on  the 
velocity  distribution  of  fragments  is  modelled  by  modifying 
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the  Gurney  equations  as  follows: 

1)  for  cylinder  charges,  continuous  casing,  it  Is 
assumed,  for  computational  purposes,  that  an 
amount  of  explosive  is  removed  from  both  ends  of 
the  cylinder.  That  amount  is  in  the  form  of  a 
cone,  which  enables  computation  of  local  C/M  to 
be  done  easily. 

ii)  for  end  projector  type  warheads,  continuous  casing, 
it  is  assumed,  for  computational  purposes,  that 
an  amount  of  explosive  is  removed  from  the  edges. 

It  is  assumed  that  only  a  cone  of  explosive  remains. 
The  local  C/M  can  then  be  computed  easily, 
iii)  for  warheads  made  of  pre-formed  fragments,  it  is 
assumed  that  the  actual  C/M  for  the  warhead  is 
reduced  to  some  effective  C/M  which  reflects  the 
greater  energy  loss. 

The  Taylor  equation  is  modified  in  accordance  with 
Ref  20.  Using  the  modified  formulae,  a  computer  program  was 
written  to  calculate  the  necessary  parameters,  velocity 
distribution  and  projection  angles,  and  plot  them  so  the 
results  could  be  compared  with  experimental  data  available. 
The  computer  program  is  contained  in  Appendix  C. 
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IV  Results  of  Application  of  Modified  Formulae 


Experimental  Data 

The  modified  equations  were  used  to  calculate  velocity 
distribution  and  projection  angles  for  end-initiated  cylin¬ 
ders  (L/D-2.0)  filled  with  various  explosives,  and  two 
different  C/M's  for  each  explosive.  Cylinders  made  of  pre¬ 
formed  fragments,  C/bfr».93  and  two  different  L/D'a  were  also 
investigated.  The  results  were  then  compared  to  the  exper¬ 
imental  data  for  those  configurations.  Experimental  data 
are  difficult  to  obtain  for  exploding  warheads.  Velocities 
can  be  measured  by  different  methods ,  such  as  X-Ray  photog- 
^  raphy,  or  electronic  screens,  but  velocities  obtained  from 

different  methods  are  slightly  different  for  the  same 
explosive/metal  geometry.  Hence,  some  methods  are  more 
accurate  than  others  and  it  is  impossible  to  give  an  exact 
value  for  the  velocity  distribution  for  a  particular  test. 
However,  it  is  assumed  that  the  data  used,  to  which  the 
equations  are  compared,  were  obtained  by  using  the  same 
measuring  techniques  and  apparatus  for  all  the  tests,  since 
they  were  obtained  by  the  same  organization.  Therefore, 
although  the  velocities  measured  may  not  be  the  exact 
velocities,  the  data  has  the  same  error  sources  due  to  the 
method  of  measurement  used.  The  data  plotted  in  Fig  11  to 
Fig  16  are  for  tests  in  which  the  geometries  were  kept  the 
(  same  for  three  different  explosives  (Octol,  Comp  B,  and  TNT) 
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and  for  two  different  C/M's  for  each  explosive.  For  Octol 
and  TNT,  the  tests  with  cylinders  of  L/D-2.0  were  repeated 
three  times  each  for  the  same  C/M  and  for  TNT,  they  were 
done  six  times  in  one  instance  and  five  in  the  other  for 
the  same  C/M.  Further  insight  can  be  gained  into  the 
accuracy  of  the  data  as  follows;  as  can  be  seen  in  Fig  11 
to  Fig  16  different  results  are  obtained  in  successive 
tests  using  the  same  C/M,  L/D  and  explosive.  For  example 
for  TNT,  L/D-2,  C/M-.77,  the  difference  is  approximately  17 
percent  between  results  of  RD899  and  RD1208;  for  Coup  B, 
C/M*. 39,  it  is  as  much  as  137.  at  some  positions  for  RD1211 
and  RD1212.  Also,  for  any  one  test,  the  data  points  do  not 
form  a  smooth  curve.  The  more  tests  of  any  particular 
configuration  (e.g.TNT,  L/D»2,  C/M-.77),  the  more  the  data 
are  spread.  The  variations  in  the  data  are  due  to  a  number 
of  factors.  First,  there  will  be  interaction  in  the  frag¬ 
ment  cloud  so  that  some  fragments  will  collide  with  each 
other  and  thus,  the  velocity  distribution  will  not  be  a 
smooth  curve  and,  depending  on  the  number  of  interactions, 
since  the  metal  casing  will  fragment  differently  in  each 
test,  there  will  be  a  difference  in  the  velocities  from 
test  to  test.  Second,  although  the  external  geometry  has 
been  kept  the  same,  the  properties  of  the  explosive  may 
differ.  The  rate  of  detonation  of  an  explosive  is  influ-  . 
enced  by  the  temperature,  pressure,  density  of  packing  and 
humidity.  Consequently,  unless  all  external  conditions  are 
the  same  for  each  test,  different  results  will  be  obtained 


and  when  all  these  differences  are  added  together,  the 
velocities  obtained  differ  between  tests  of  same  L/D,  C/M, 
and  explosive.  Although  the  velocity  distribution  has  a 
certain  amount  of  error,  most  of  the  tests  (except  for  some 
tests  with  TNT)  show  the  influence  of  end  effects. 

Cylinders  with  Continuous  Casings 

Figures  11  to  16  show  the  results  from  the  Gurney, 
modified  Gurney  and  modified  Taylor  equations  plotted  with 
the  experimental  data  for  explosive  cylinder  tests,  for 
each  case,  L/D-2.0.  Because  the  experimental  data  are 
different  for  tests  of  the  same  explosive,  L/D  and  C/M,  the 
estimates  from  the  Gurney,  modified  Gurney,  and  modified 
Taylor  equations  should  be  compared  to  an  average  of  the 
test  results  for  each  explosive  and  C/M  for  a  more  meaning¬ 
ful  comparison. 

The  upper  graph  of  Fig  11  shows  the  velocity  distribu¬ 
tion  obtained  by  the  Gurney  equation  (dashed  line)  and  the 
modified  Gurney  equation  (solid  line)  plotted  with  the  exper¬ 
imental  data  for  cylinders  filled  with  Octol,  C/M-.86.  It 
can  be  observed  that  the  Gurney  equation  estimate  reasonably 
well  fragment  velocities  between  .45Sx/L*.85.  However,  the 
estimate  is  progressively  larger  in  error  as  .452x/LiO  and 
.85ix/LS1.0.  On  the  other  hand,  the  estimates  from  the 

modified  Gurney  equation  are  reasonably  accurate  for  •  15 -jl-  .  95k 

l 

The  modified  Gurney  equation  is  reasonably  accurate  for 
approximately  80  percent  of  the  relative  distance,  whereas. 
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the  Gurney  equation  is  only  accurate  for  40  percent  of  the 
relative  distance.  The  modified  equation  reproduces  fairly 
closely  tests  RD875  and  RD1217  between  .  lix/14.95. 

The  bottom  graph  compares  the  results  of  the  modified 
Taylor  equation  vith  the  experimental  data  for  the  same 
cylinder  tests.  As  can  be  observed,  the  modified  Taylor 
equation  gives  a  reasonably  accurate  estimate  compared  to 
the  estimate  from  the  normal  Taylor  equation  (see  Fig  7) . 

The  superiority  of  the  modified  Taylor  equation  over  the 
normal  Taylor  equation  is  clearly  evident  in  this  case. 

In  the  upper  graph  of  Fig  12,  the  velocity  distribution 
obtained  by  the  Gurney  and  modified  Gurney  equations  are 
plotted  with  the  experimental  data  for  cylinders  filled  with 
Comp  B,  C/M*«.79.  It  can  be  observed  that  the  modified 
Gurney  equation  gives  a  reasonably  accurate  estimate  of  the 
velocity  distribution  for  . lix/L^.95,  whereas  the  Gurney 
equation  is  only  accurate  for  .45x/Lf.85.  Furthermore  the 
modified  equation  reproduces  fairly  closely  test  RD588  be¬ 
tween  . l*x/L*.95. 

The  bottom  graph  compares  the  results  of  the  modified 
Taylor  equation  with  the  experimental  data  for  the  same 
cylinder  tests.  As  can  be  observed,  the  modified  Taylor 
equation  gives  a  reasonably  accurate  estimate  compared  to 
the  estimate  from  the  normal  Taylor  equation  (see  Fig  7) . 

The  superiority  of  the  modified  Taylor  equation  over  the 
normal  Taylor  equation  is  clearly  evident  in  this  case. 
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Th«  im  infonution  is  contained  in  Fig  13  for 
cylinders  filled  with  TNT,  C/M-.77.  In  both  graphs,  the 
data  show  a  wide  spread  and  the  agreement  of  the  equations 
with  the  data  is  not  as  good  as  in  the  previous  cases. 
However,  the  modified  equation  reproduces  fairly  well  test 
RD1208.  Most  of  the  tests  show  that  end  effects  axe  present 
and,  for  this  reason,  the  modified  Gurney  equation  is  as¬ 
sumed  to  give  better  estimates  than  the  Gurney  equation. 

The  estimates  of  the  projected  angles  are  clearly  too  high. 
However,  the  modified  Taylor  equation  is  superior  to  the 
normal  Taylor  equation. 

Figure  14  shows  the  velocity  distribution  obtained 
from  the  modified  Gurney  and  Gurney  equations  in  the  upper 
graph  for  cylinders  filled  with  Octol,  C/M-.43.  Again,  it 
can  be  observed  that  the  modified  Gurney  equation  gives 
better  estimates  than  the  normal  Gurney  equation.  The  bottom 
graph  shows  the  estimates  of  the  projected  angles  from  the 
modified  Taylor  equation  and  experimental  data  for  the  same 
tests.  The  agreement  of  the  modified  Taylor  equation  with 
the  experimental  data  is  better  than  that  with  the  normal 
Taylor  equation. 

Figures  IS  and  16  contain  the  same  information  as  the 
other  figures  but  for  cylinders  filled  with  Comp  B,  C/M-.39, 
and  cylinders  filled  with  TNT,  C/M-.38,  respectively.  In  Fig 
IS,  agreement  of  the  equation  with  the  experimental  data  is 
not  as  good  as  the  agreement  in  most  of  the  previous  figures. 
However,  the  modified  equations  are  superior  to  the  normal 
equations  for  obtaining  estimates  of  velocity  and  projected 


angle .  The  experimental  data  in  Fig  16  shows  a  very  large 
spread.  Both  the  Gurney  and  the  modified  Gurney  equations 
overestimate  the  velocity  distributions.  Neither  equation 
gives  the  good  agreement  obtained  in  most  of  the  other  cases. 
The  modified  Gurney  equation  should  be  superior  to  the 
normal  Gurney  equation  because  it  accounts  for  end  effects, 
which  although  not  evident  from  the  data,  are  present. 

The  results  of  the  modified  Gurney  equation  are  not 
plotted  for  x/L<.15  and  x/L>.95  because  the  velocity  de¬ 
creases  rapidly  to  zero  in  those  regions  and  this  contradicts 
reality  since  the  end  fragments  will  have  a  finite  velocity. 

A  better  representation  would  be  to  continue  the  fragment 
velocity  curve  with  the  same  constant  slope  obtained  for 
the  last  points  in  the  plot.  However,  since  there  are  few 
experimental  data  ppints  in  those  regions,  it  is  difficult 
to  ascertain  whether  this  procedure  should  be  followed. 

Since  the  modified  equation  gives  reasonable  agreement  in 
the  region  . 15<x/L<. 95 ,  it  is  deemed  appropriate  to  repre¬ 
sent  end  effects . 

Cylinders  with  Pre- Formed  Fragments 

The  equations  are  valid  for  cylinders  made  with  pre¬ 
formed  fragments.  However,  since  the  fragments  separate 
earlier  than  cylinders  with  continuous  casings  and  thus  let 
the  energy  dissipate  faster,  it  is  expected  that  the  veloci¬ 
ty  should  be  lower.  This  is  confirmed  by  the  experimental 
data  plotted  in  Fig  17.  Although  the  C/M  is  .93,  the  maxi¬ 
mum  fragment  velocity  attained  is  approximately  seven 
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percent  lover  than  tests  using  the  same  explosive  and  L/D. 
but  where  the  C/H  is  only  .86  (see  Fig  11).  If  no  modifi*- 
cation  for  pre-formed  fragments  is  made,  the  results  of  the 
Gurney  equation  are  in  error  by  more  than  10  percent,  as 
shown  in  the  upper  graph  of  Fig  17.  The  lower  graph  shows 
the  results  of  the  modified  Taylor  equation.  Since  there  is 
less  energy  to  accelerate  the  fragments,  this  is  equivalent 
to  using  a  smaller  effective  C/M.  Consequently,  starting 
with  the  results  of  the  Gurney  equation,  modified  for  end 
effects ,  the  C/M  of  the  charge  was  reduced  until  good  agree¬ 
ment  with  experimental  data  was  obtained.  This  occurs  when 
the  effective  C/M  is  taken  to  be  80  percent  of  the  actual 
C/M.  Using  this  effective  value  for  C/M  in  the  modified 
Gurney  equation  yields  the  results  shown  in  the  upper  graph 
of  Fig  18.  The  C/M  used  in  the  normal  Gurney  equation  has 
also  been  reduced  by  20  percent.  The  difference  between 
the  modified  equation  and  experimental  data  is  on  the  order 
of  one  percent  for  the  velocity  distribution  for  .l<x/L<.95. 
The  plot  in  the  bottom  graph  gives  a  fairly  good  represen¬ 
tation  of  the  projected  angles.  Tables  VII  and  VIII  give 
the  numerical  results  of  the  standard  and  modified  equations. 

When  cylinders  with  L/D<2  are  investigated,  the  modi¬ 
fied  equations  give  slightly  less  satisfactory  results. 
Experimental  data  for  tests  with  cylinders  of  L/D-1.0  are 
compared  to  the  results  of  the  modified  equations  in  Fig  19 . 
The  same  corrections  as  those  for  cylinders  with  L/D-2.0 
are  used.  As  can  be  seen,  the  equation  slightly  overesti¬ 
mates  the  velocities  at  x/Le.7  and  underestimates  them 
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towards  Che  ends.  The  maximum  fragment  velocities  attained 
are  more  than  20  percent  lower  than  that  for  L/D-2.0,  even 
though  the  C/M  is  .93  for  both  cases.  This  is  because  the 
rarefaction  waves  have  a  greater  effect  on  the  cylinder  and 
cause  a  greater  loss  in  velocity.  Using  the  same  modifications 
as  that  for  cylinders  with  L/D*2.0  give  a  larger  error  (approx¬ 
imately  9  percent  at  x/Le.2)  but  the  equation  gives  acceptable 
results  for  .  2<x/L*.9  which  are  clearly  superior  to  the  esti¬ 
mates  given  by  the  normal  Gurney  equation  corrected  for  pre¬ 
formed  fragments  only.  A  different  correction  could  be  used 
to  obtain  better  results,  but,  since  the  results  are  fairly 
good,  the  same  correction  factors  should  be  used  for  all 
cylinders  of  different  L/D's  to  retain  a  simple,  easy  to 
use  method.  Table  IX  gives  the  numerical  results  of  the 
standard  and  modified  equations. 

From  the  results  obtained,  it  can  be  seen  that  the 
Gurney  equation  for  cylinders  made  with  continuous  casings , 
modified  for  end  effects  as  shown  in  Fig  9,  will  give 
reasonably  accurate  estimates  of  the  fragment  velocity 
distribution.  Also,  for  cylinders  made  with  pre- formed 
fragments,  the  results  show  that  the  C/M  should  be  re¬ 
duced  by  20  percent,  in  addition  to  the  modification  for 
end  effects. 

End  Projectors 

The  results  of  the  corrected  Gurney  equation  for  end 
projectors  are  tabulated  in  Table  X.  Since  pre-fomied 
fragments  are  used,  the  effective  C/M  used  in  making  the 
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FIG  11.  VELOCITY  DISTRIBUTION  RND  PROJECTED  RNGLE 
FOR  CYLINDER  FILLED  MITH  OCTOL.  L/D=2.0,  C/f1=.86 
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FIG  15. VELOCITY  DISTRIBUTION  RND  PROJECTED  ANGLE 
FOR  CYLINDER  FILLED  NITH  COMP  B.  L/0=2.0.  C/H=.39 


Table  VI 


r 


ot«oeomooMNvooN<9ni-iintQoinvfi 

r-r>>cnoMnr^vo«oooooo»Ai-tvovO'i-eso4*-t 

NOHHNM«(MnoHM<n(niN®«nmo 

irimtM  i-i(Nfnn>jirufiininimrnoeooN 


C  CJ 

c 

<U 

8 

o 

5>sO 

* 

o 

«t 

>nH 

N 

vss 

►» 

6s 

0) 

r* 

3,C 

fi 

O  u 

3 

H 

•rl 

O 

01 

•o  3 

> 

cd*0 

u  oo 

•t)  01 

a 

Sp* 

C> 

Op-1 

OB 

M  S 

■Ut» 

3. 

g 

cob 

A 

T)  M 

C  V 

1 

s 

3*o 

«n 

_  c 

ON 

•o^ 

• 

01  r~l 

VO 

i-H 

■H  >N 

1 

0)  0) 

*W  U 

>N 

>  3. 

*r« 

w 

'-x. 

•o  M 

T-l 

4->  g 

5  o 

o 

t4  B 

o 

c 

rH 

M 

Hr»ver*inNo\n>fHN«j'nr^ooo\oo<n 

NrtrvNioin^'Jmff\ooinr4«)0'JvOHCo 

lAMHOMrio^fr'ff'Ot-iNnfOvincirx 

r<cNNrt^><t^'Jinininimnininin«4-rt 


^nvovowommmwoMJ'O'OMMnovoM 

NcMHr(oOMioNNrsNrsr^r>.inoonN 

vD»opiNNnMn«mncirtio«NHco 


ON  On  On  On  On  O'  O'  On  O'  O'  O'  O'  O'  on  on  On  On  On  O' 
fsfsNrN*  r»  n  fN»K  f''  r-  r^*  n  K  r%*  r** 

nMtnnnnwi'inrmcnorNnnFK'ifO 


ooooooooooooooooooo 

motflomomomomoi/ioinoinoin 

OHHNNronvtvfmiAvOior^rvooiacM^ 


53 


PROJ  RNG  t  DEO)  FRflO  VEL  ( rtfl/A  8 ) 


REL  01 ST  RLONO  CASINO 


8 


HO  16 . VELOCITT  DISTRIBUTION  RNO  PROJECTED  RNOLE 
FOR  CYLINDER  FILLED  WITH  THf .  L/D=2.0,  C/H=.3B 


Table  VII 


1 

4 

i 


meiFiEO 


r\ 


^S^ff'csjr^ir'fnr^CT'CMoo^<»(^-<tvocsi 

00tSO\N0\r-IHOiU1ONfiriv£lNr0vl-NN 


f"  I/O  r-H 


■dC'0'4-<t<i">lOlO<©vOvCf‘-OOi 


•a 

ss. 

O  4-1 

^  c 
a;  -c  <u 

gig, 

ls2 

•o  u, 

X)  4) 
d  rH  XJ 

Q  r-4  0) 

•O-H  8 

gh  C 

CD  O 
*J  Mh 

W  CD  | 


S2^I^^2i,3!2i'>‘0,-,ir'0',^-'-<0'vOirnN 
fiJP'^C^Cl,ri'3'<r*cs|'^'0000\OCOaOVDt^)oO<t' 
vDvovoinnai'j’ooi—i  co  ^  in  \o  r*-  n  \o  nvo>j 

•HfN4cnsf-<fLni/‘>voic>»oiOvov£>iovoioin'<a- 


XJ  0)  6 

•a  C 

g 

C  *d  pg 

'3- 

<0  i-g 

r- 

o 

-  ^  u 

. 

eg 

DU  O 

N 

a>  «w 

co 

oo 

•d  Vt 

cr> 

H 

*4-1  O  X) 

to 

■d  <4-1  Q) 

* 

4-1 

XJ  4J 

oo 

•d 

O  ° 

B 

o 

X  <U 

£ 

o 

>g 

rd 

<4-1  U 

o 

<U 

O  O 

> 

u 

c 

c 

o 

o 

a 

•d 

•d 

o 

4J 

d 

(0 

to 

§* 

CM 

A 

c 

o 

4J 

o 

0* 

a 

hJ 

a 

SSS^iS^i2?'J0riOOCS,eslCslCs,0OOOMOO 

SSJSS5^2Ji£!ll';,‘inv0v0v<:)'0,r>,“,o'ONn4 

OvcnnvOCOOOO'OOOOOOOOOOOvD?! 

»^M^^r-li--IC'l«SlevlC\llSltNlC4C\lr>4r-l\--Ir-4 


Sf^f)!f^eiJCi!tNCV,'NCN,CN,ev,CNICSICVI£^nlcs4CM 

^^^^^D^OiOvOvOvOvOvOiOiOiOiDvOiDiO 

ooooooooooooooooooo 

<n  CNJ  (N|  cvj  C-4  CVJ  es  CN  cs  csi  C\|  cvi  cvi  C\1  eg  eg  <N  C>j  cvi 


OOOOOOOOOOOOOOOOOOO 

jo^oinoinoiooinoinoifioinom 

Or-trgfgrgcncn-s-^imnvoioSrCSSSS 


18.864 


>\  to 

QJ  0) 
*“•£3  i 


3-»~«  '  ~~ZZ2°"'e-r*8lS 2 

1  *  i  i  i 


„*  *  •  •  ■ ’ 


Sa„ 

g  °£ 

|*S 

i  "gi 

I  S35 

1  I?? 

Sfc§ 

«  fcfc 
W«Y 

"2  «i  , 

§3*  ? 

i 

zt-v  l 
« 


60 

4f  60 

r“i 


»  OJ 

00  : 

Ox  OQ 

I'  V 

<o  i? 


ss|g52l5ss|aga|2^sg 


«55»?ssi535Ssta3SSo„ 

^^“ZSXSig 


isSiiissgggaa, 


o  i, 

<8 

^  C 

5  I 
5  * 


°^mm§  sgggggg 
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calculations  was  reduced  by  20  percent.  Figure  20  shows 
the  geometry  of  the  end  projector  with  a  comparison  of  the 
results  obtained  from  the  Gurney  equation  (Eq  22)  and 
experimental  data  (Ref  9).  If  no  corrections  are  made,  the 
Gurney  equation  gives  a  maximum  velocity  of  approximately 
2.7  mm/ps.  The  experimental  data,  although  it  may  not  give 
the  exact  velocity  for  reasons  discussed  previously,  give 
a  maximum  velocity  of  approximately  2.1  mm/ps,  a  difference 
of  22  percent  from  that  of  Eq  22.  However,  when  the  modified 
Gurney  equation  (Eq  32)  is  used,  the  difference  is  only  11 
percent.  Unfortunately,  there  is  only  one  set  of  experimental 
data  to  test  the  modified  equation.  With  more  experimental 
data,  it  might  be  possible  to  derive  a  simple,  more  precise 
modification  factor.  Nevertheless,  the  present  modification 
is  simple  and  gives  reasonable  results.  It  will  give 
results  of  the  same  form  for  any  other  L/D  and  C/M,  since 
changing  any  one  or  the  other  ratio  changes  the  size  of  the 
cone  of  the  explosive  accordingly.  The  modified  equation 
should  therefore  give  results  which  are  comparable  to  the 
experimental  data. 

A  slightly  different  modification  factor  for  the 
Gurney  equation  for  an  open-faced  sandwich  was  found  in 
Ref  14*  The  factor  proposed  is  the  removal  of  a  cone  of 
explosives  with  an  angle  of  30  degrees  from  the  normal  at 
the  edges  of  the  charge  (see  Fig  20). 

The  results  obtained  if  that  factor  is  used  are 
shown  as  the  dotted  line  in  Fig  21  and  in  Table  X. 
Unfortunately,  only  one  set  of  experimental  data  was  found 


Fig  20.  Correction  for  End  Eff erf.  «_ 

rfects  for  Open-Faced  Sandwi< 

“  C7a"  tHe  0,0  A‘  -  b.  seen,  the  results 

“  ne<i  "*  nCt  ^  80°d  “  «—  Siven  by  the  Bdl£iMtta 
factor  given  la  „  32.  On  the  basis  o£  that  eitparlmantal 

data,  the  -edification  proposed  in  Section  III  of  this 
study  is  superior  to  the  one  proposed  in  Ref  14. 

In  general,  given  a  cylindrical  warhead  or  an  end 

projector  type  warhead  of  a  certain  L/D  and  C/M.  a  reasonab! 

accurate  estimate  of  the  velocity  distribution  for  the 

fragments  can  be  obtained  by  us  in?  r 

y  g  the  Gurney  equations  with 

appropriate  modifications  for  the  narticl 

cne  particular  geometry  in¬ 
volved. 

The  angle  of  departure  of  the  fragments  was  not  found 
acauee  it  would  be  a  very  compler  undertaking.  The  Taylor 
ormu  a  applies  to  a  system  where  the  detonation  wave  is 
perpendicular  to  the  surface  to  be  projected,  if  the  charge 
in  an  end  projector  is  point-initiated,  a  spherical  detons. 
«°n  wave  will  be  formed.  Reflection,  will  occur  at  free 
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surfaces  and  Che  pattern  of  Che  detonation  wave  will  be  quite 
difficult  to  obtain.  An  analysis  similar  to  Taylor's  could 
be  done  for  this  situation  but  since  the  tocal  dispersion  of 
the  fragments  is  of  greater  interest  than  that  of  individual 
fragments,  no  attempt  was  made  to  obtain  the  projection  angle 
of  the  fragments.  The  total  dispersion  can  be  controlled  by 
imbeding  the  fragments  in  a  plastic  ring  as  was  done  in 
the  experiments  (Ref  2)  which  will  be  discussed  next. 

An  attempt  was  made  to  use  the  method  on  a  different 
type  of  end  projector  (see  Fig  22).  As  can  be  seen,  the 
projector  consists  of  five  layers.  The  uppermost  is  comprised 
of  32  aluminum  cubes  surrounded  by  plastic.  The  other 
layers  are  C-4,  steel,  1/16  inch  thick  Detasheet,  and  plywood, 
in  that  order.  The  projector  has  a  12.7  cm  diameter.  The 
Gurney  equation  for  an  asymmetrical  sandwich  was  used  for 
this  analysis.  Experimental  results  are  contained  in  Ref  2. 

As  is  explained  in  the  reference,  measurement  of  the  fragment 
velocities  was  not  easy.  Different  methods  were  used  and, 
then,  only  the  speed  of  the  fastest  fragment  could  be 
obtained  for  reasons  outlined  in  Ref  2.  Although  the  tests 
were  with  essentially  the  same  L/D's,  C/M's,  and  explosive 
fillings,  the  measured  velocities  differed  greatly. 
Consequently,  calculated  results  are  best  compared  to  the 
average  measured  velocity  for  any  one  series  of  tests. 

Since  only  the  fastest  fragment  speed  was  recorded,  which 
should  be  the  one  in  the  middle  of  the  circular  disk,  the 
results  of  the  Gurney  equation,  with  only  a  correction  for 
using  pre-formed  fragments,  are  compared  to  the  experimental 
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data.  Detailed  calculations  for  one  specific  test  (Shot 
No  43)  are  shown  in  Appendix  B. 

The  calculated  velocity  from  Gurney’s  equation  for  Shot 
43  is  2200  m/s,  and  the  average  measured  speed  of  the  fastest 
fragment  for  that  series  of  tests  is  2430  m/s.  Thus,  the 
calculated  value  underestimates  the  velocity  by  approxi¬ 
mately  ten  percent.  The  calculations  were  repeated  for 
Shot  No  66  in  which  almost  twice  as  much  C-4  was  used. 

The  velocity  was  calculated  to  be  2690  m/s,  and  the 
average  measured  velocity  for  that  series  of  tests  is 
3410  m/s.  In  this  instance,  the  method  underestimates  the 
velocity  by  about  21  percent.  At  this  point,  it  seems  that 
the  modified  method  should  not  be  applied. 

The  method  does  not  work  well  because  these  end  pro¬ 
jectors  are  not  end  initiated.  The  method  of  initiation  is 
important  in  what  occurs  during  the  detonation  process. 

Since  the  C-4  is  initiated  simultaneously  around  the  peri¬ 
meter  by  the  Detasheet,  a  circular  shock  wave  will  form  and 
converge  towards  the  center  of  the  charge.  The  shock  wave 
interacts  with  itself  to  produce  an  amplified  wave.  The 
results  are  similar  to  a  cylinder  which  is  initiated  at 
both  ends.  Such  a  case  is  considered  in  Ref  18.  The  velocity 
of  the  fragments  at  the  middle  is  approximately  15  percent 
greater  in  such  instances  than  when  the  cylinder  is  initiated 
from  one  end  only.  This  implies  that  no  reduction  should  be 
made  in  the  C/M  because  of  pre-formed  fragments,  or  that  the 
modification  can  be  made  and  the  results  increased  by  15 
percent.  Since  no  additional  information  is  known,  the  value 
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1*  approximate  because  different  configurations  might  yield 
different  percent  increases.  However,  this  value  should  be 
close  to  any  value  obtained  for  a  different  configuration. 
Consequently,  if  the  total  effective  C/M's  are  not  reduced 
in  Gurney's  formula,  the  results  obtained  are  velocities  of 
2375  m/s  and  2870  m/s  for  Shots  43  and  66  respectively. 

This  reduces  the  error  to  approximately  two  percent  for  Shot 
43  and  16  percent  for  Shot  66.  It  can  be  concluded  that, 
when  initiation  causes  amplification  of  the  detonation  wave, 
such  as  simultaneous  initiation  of  a  cylinder  from  both  ends 
or  initiation  as  in  this  type  of  projector,  no  correction 
should  be  made  for  using  pre-formed  fragments.  It  may  be 
that  the  C/M  should  be  increased  a  certain  amount,  but 
insufficient  data  precludes  that  conclusion. 

Once  again,  obtaining  the  deflection  angle  of  the  frag¬ 
ments  is  a  very  difficult  task.  From  the  experiments  presen¬ 
ted  in  Ref  2,  it  is  evident  that  having  the  plastic  ring  (see 
Fig  22)  around  the  cubes  cut  down  dispersion  greatly.  The 
plaatic  ring  had  a  12.7  cm  diameter  and  the  cubes  were  5.72  cm 
in  length.  In  these  instances,  the  solid  angle  for  75  percent 
of  the  fragments  at  a  distance  of  24  feet  ranged  from  1.3  to 
9.1  millisteradians .  It  appears  that,  in  general,  by  enclosing 
the  intended  fragments  In  a  similar  plastic  ring  whose  diameter 
is  twice  the  length  of  the  fragments  arranged  in  a  square, 
similar  dispersions  should  be  obtained.  Therefore,  in  design¬ 
ing  an  end  projector  type  warhead,  the  fragments  should  be  en¬ 
closed  by  some  sort  of  ring,  if  a  small  dispersion  is  desired. 
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32.  A1  Cubes 


V  Conceptual  Warhead  Designs 


Analysis  of  conceptual  preliminary  designs  of  war¬ 
heads  can  now  be  made  using  the  information  gained  in 
the  previous  sections.  For  example  purposes,  warhead 
designs  which  will  give  initial  fragment  velocities  of 
3000  m/s  will  be  considered.  Also,  it  is  assumed  that 
the  warhead  is  to  be  fitted  in  an  air-to-air  missile 
and  is  detonated  at  a  distance  of  15  meters  from  the  target. 

Velocity  of  Fragments  at  Target 

It  is  important  to  know  what  the  final  velocity  of  the 
fragments  will  be  when  they  reach  the  target.  If  they  are 
too  low,  that  is,  if  the  kinetic  energy  is  low,  they  will 
not  penetrate  the  target.  Since  the  warhead  is  not  likely 
to  be  fitted  at  the  nose  of  the  missile,  the  fragments  will 
have  to  go  through  some  of  the  missile's  components  when 
the  warhead  detonates.  Consequently,  the  fragments  will 
leave  the  missile  at  a  lower  velocity  than  3000  m/s ,  not 
Including  the  missile's  own  velocity  which  is  imparted 
to  the  fragments.  The  warhead  is  made  of  pre-formed  steel 
cubes  (.953  cm)  and  it  is  assumed  that  they  would  go  through 
a  number  of  components  which  are  comparable  to  a  2.54  cm 
al  u minim  sheet  (exact  thickness  would  be  known  in  an  actual 
design).  If  plugging  type  failures  are  assumed,  the  residual 
velocity  of  the  fragments  can  be  found. 
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-  - ..  A-  r  v2 

R  Mp  +  M0  L  VP 


(Ref  1) 


(33) 


where 

VR  -  residual  velocity  (m/s) 

Vp  -  initial  fragment  velocity  (m/s) 

Mp  *  mass  of  steel  fragment  (-6.7  g) 

-  mass  of  Al  plug  punched  by  the  fragment 
^  (-6.4  g) 

V.  -  ballistic  limit,  velocity  at  which  50%  complete 
penetration  and  50%  partial  penetration  of  the 
target  plate  can  be  expected 

An  approximate  numerical  value  for  the  V5Q,  ballistic  limit, 
can  be  found  by  using  an  equation  obtained  to  fit  . 30  cal 
armor-piercing  projectiles  because  an  exact  value  could  not 
be  found  from  the  sources.  Assuming  that  the  fragments  go 
through  2024-T3  Al, 


A 


thickness  of  plug 
plug  diameter 


-  575.5  m/s 


(Ref  3)  (34) 


VR  "  T371  (3000*  -  575. 52  )** 


1506  m/s 


(35) 


The  velocity  at  15  m  from  the  point  of  detonation  can 
be  found  by  using  an  aerodynamic  drag  law 

CnpA 


Vn  exp 


(Ref  1) 


(36) 


where 
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VM-  fragment  velocity  (m/s) 

x  -  distance  from  point  of  detonation  (-15  m) 

Vq  -  initial  velocity  (-1506  m/s) 

Cp  -  drag  coefficient  (-1.6b) 
p  -  air  density  (-.00123  g/cc) 

j 

A  -  fragment  presented  area  (-.908  cm  ) 

M  -  fragment  mass  (-6.7  g) 

For  rotating  cubes,  was  obtained  from  Ref  7.  For  compu¬ 
tational  purposes,  it  is  assumed  that  the  warhead  is  detonated 
at  sea  level.  Therefore,  the  fragments  will  have  a  velocity 
of 

VM  -  1224  m/s  (37) 

However,  due  to  end  effects,  the  fragments  towards  the  edges 
will  have  lesser  velocities.  From  experimental  data  for  the 
first  end  projector  discussed  in  the  previous  section,  the 
velocity  of  the  fragments  located  at  a  distance  of  x/L-.25 
from  the  edge  i3  reduced  by  approximately  17  percent.  In 
this  case,  these  fragments  would  attain  a  velocity  of  approxi¬ 
mately  1030  m/s. 

These  velocities,  given  favorable  relative  velocities 
between  the  target  and  the  missile,  should  be  sufficiently 
high  to  penetrate  thin  aircraft  skins  and  components  not 
protected  by  armor. 

If  a  sufficient  number  of  fragments  pierce  and  destroy 
a  vulnerable  area  (cockpit,  engines,  fuel  tanks,  etc...),  the 
target  will  be  destroyed.  Therefore,  by  using  the  modified 
Gurney  equation  for  an  end  projector  type  warhead  (Eq  32) , 
one  finds  that  the  velocity  of  approximately  50  percent  of 
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the  fragments  will  be  greater  than  1030  m/s,  a  velocity 
deemed  sufficiently  high  to  destroy  a  target.  However, 
sources  were  not  found  to  confirm  this.  Experimental  testing 
would  have  to  be  carried  out  for  fragments  at  that  velocity 
against  targets  of  interest.  For  illustrative  purposes,  it 
is  assumed  that  a  warhead  which  imparts  an  initial  velocity 
of  3000  m/s  to  the  fragments  near  the  middle  of  the  warhead 
surface  would  be  desirable.  Conceptual  designs  of  different 
geometries  which  produce  fragments  at  3000  m/a  will  be 
investigated. 

Example  1 

The  first  design  of  an  end  projector  is  a  simple  one 
using  Comp  B  as  the  explosive  charge  and  32,  .953  cm  steel 
cubes  enclosed  in  a  plastic  ring  in  a  configuration  like  the 
one  shown  in  Fig  10  and  Fig  20.  Based  on  similar  projector 
tests  seen  in  the  previous  section,  enclosing  fragments  in 
a  plastic  ring  should  give  a  small  dispersion  at  15  m.  Also, 
the  warhead  is  given  a  12.7  cm  diameter  to  fit  inside  an 
air-to-air  missile. 

To  find  the  characteristics  of  the  warhead,  Gurney's 
equation  for  an  open-faced  sandwich  is  used  (Eq  22).  Solving 
for  C/M,  with  «  3000  m/s  ,  the  warhead  will  have  a 
C/M  ■  4.47  .  The  other  root  of  the  equation  is  negative  and 

is  therefore  neglected.  Since  pre-formed  fragments  are  used, 
the  C/M  must  be  increased  by  20  percent,  or  to  C/M-5.36  , 

to  compensate  for  the  loss  of  energy  between  the  separating 
fragments.  Also,  it  can  be  observed  that,  from  the  first 
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end  projector  tests,  the  method  still  overestimates  the 
velocity.  To  reproduce  the  data  points  at  the  center  of 
the  graph  (Fig  20)  more  exactly,  a  value  of  C/H  -  1.85  has 
to  be  used  for  that  experiment  rather  than  C/M  “2.94  , 

the  actual  value.  This  corresponds  to  a  decrease  of  approx¬ 
imately  37  percent  of  the  actual  C/M.  Consequently,  to  obtain 
the  most  exact  estimate  of  the  velocity  in  this  first  design, 
the  C/M  to  be  used  in  the  equation  should  be  C/M  -  4.47  x 
1.37  -  6.12  . 

The  total  weight  of  the  steel  cubes  is  M  -  214.4  g 

9 

The  weight  of  the  plastic  ring  (same  composition  as  the  one 
used  in  Ref  2)  is  Mp  “  156  g 

Total  weight  of  upper  level,  M  -  Mg  +  Mp  “370. 4g 
Therefore,  the  required  weight  of  Comp  B  is  C  “  2266.9  g  . 

The  density  of  Coup  B  is  p  -  1.75  g/cc  (Ref  4). 

The  height  of  the  Comp  B  cylinder  is  10.2  cm. 

Therefore,  a  warhead  with  L/D  -  .881  and  C/M  -  6.12  would 
have  a  velocity  distribution  with  the  center  fragments  at 
approximately  3000  m/s  (see  Fig  24a).  The  C/M  for  this  warhead 
is  quite  high.  It  even  falls  outside  the  range  of  applica¬ 
bility  of  the  Gurney  theory.  This  warhead  would  not  be  selected 
because  its  high  C/M  would  probably  cause  the  fragments  to 
break  up. 

Example  2 

A  slightly  more  complex  example  is  presented.  The 
effects  of  the  addition  of  a  tamper  weight,  or  driver  plate, 
to  the  first  design  are  investigated.  The  configuration  of 
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the  warhead  is  the  same  as  that  of  an  asymmetric  sandwich 
and  Gurney's  equation  (Eq  20)  for  that  geometry  is  used  to 
find  the  dimensions  of  the  components  of  the  warhead,  L/D, 
and  C/M. 

Although  there  does  not  appear  to  be  a  physical  optimum 
for  C/M  and  N/M  for  a  particular  design,  it  is  possible  to 
use  the  equation  to  obtain  a  "best"  design  for  a  given  set 
of  specifications,  such  as  fragment  velocities  to  be  attained, 
maximum  warhead  size,  and  production  cost  of  the  warhead, 
etc...  The  results  of  the  equation  are  tabulated  in  Table  XI 
for  a  few  values  of  C/M  and  N/M.  The  non-dimensional  velocity 
can  be  obtained  from  it,  and  given  an  explosive,  i,e,  /TE~  , 
the  fragment  velocity  can  be  obtained  for  that  particular 
C/M  and  N/M. 

From  the  table ,  it  can  be  seen  that  as  C/M  and  N/M 
increase,  the  velocity  increases.  Although  the  largest 
N/M  used  was  10,  larger  values  give  increasing  velocities. 

The  table  was  assembled  for  .liC/MilO  ,  but  since  the 
Gurney  equations  give  best  results  for  ,1<C/M<5  ,  the 

velocities  obtained  for  C/M  -  10  are  largely  approximate 
and  are  tabulated  only  to  show  that  V  /ZF  does  in  fact 
increase  as  N/M  and  C/M  increase.  The  data  are  plotted  in 
Fig  23.  The  solid  lines  indicate  increasing  velocities 
with  increasing  C/M  and  N/M.  The  dotted  lines  represent 
a  warhead  where  the  total  weight  was  kept  constant  at  3  kg. 

N  and  M  were  varied  while  C  was  kept  constant;  this  was 
repeated  for  three  different  values  of  C.  It  can  be  observed 
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Table  XI 

Values  of  V^/ZSIT  for  Various  C/M  and  N/M 


\C/M 
N/M  X. 

0.1 

0.5 

1.0 

2.0 

5.0 

10.0 

0.1 

0.117 

0.355 

0.558 

0.818 

1.175 

1.392 

0.5 

0.185 

0.427 

0.608 

0.837 

1.167 

1.379 

1.0 

0.222 

0.480 

0.655 

0.866 

1.168 

1.369 

3.0 

0.269 

0.536 

0.745 

0.943 

1.196 

1.361 

5.0 

0.284 

0.593 

0.781 

0.982 

1.223 

1.370 

10.0 

0.311 

0.655 

0.866 

1.095 

1.270 

1.400 

that  for  a  fixed  total  weight  and  a  fixed  amount  of  explosive, 
the  velocity  increases  as  the  tamper  weight,  N,  increases 
and  the  fragment  mass,  M,  decreases.  Also,  the  larger  the 
amount  of  explosive,  the  higher  the  velocity.  These  curves 
would  be  different  for  a  different  total  weight.  They  are 
useful  in  designing  a  warhead.  For  instance,  if  the  speci¬ 
fications  called  for  a  warhead  with  a  maximum  weight  of  3  kg, 
fragment  velocities  of  at  least  2000  m/s,  and  Comp  B  as  the 
explosive  to  be  used,  the  graph  would  give  the  different 
combinations  possible  which  would  satisfy  those  specifications. 
In  this  case,  those  combinations  are: 


C  (kg) 

C/M  . 

N/M 

M  (kg) 

N  (kg) 

0.5 

@  1 

@  4 

0.5 

2.0 

1.0 

@  1.25 

@  1.5 

0.8 

1.2 

1.5 

@  1.5 

@  0.6 

1.0 

0.5 
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The  selection  could  then  be  made  on  which  of  these 
combinations  meets  best  the  specifications. 

Returning  to  the  second  example  where  a  fragment 
velocity  of  3000  m/s  is  required  and  Comp  B  is  to  be  used, 
the  solid  lines  in  Fig  23  can  be  used  to  find  possible  combi¬ 
nations  to  meet  the  velocity  of  V^/ £  1 . 12 ,  i.e.  V^/ZZE” 
equals  3000/2680  *  1.12  .  Some  combinations  are  : 


C/M 

N/M 

2.5 

20.0 

3.0 

7.0 

3.5 

4.0 

4.0 

2.5 

More  possibilities  exist  but  for  C/M<2.5  ,  large  N/M's 

are  required,  which  would  probably  make  the  design  more  diffi¬ 
cult,  and  for  C/M>4  ,  the  Gurney  equation  starts  to  lose 

its  applicability.  The  values  of  C/M  «  3.0  and  C/M  -  4.0 
will  be  used  design  the  warhead;  they  can  be  compared  to 
find  which  design  is  likely  to  be  the  best.  Since  pre¬ 
formed  fragments  are  to  be  used,  the  C/M  ratio  must  be 
increased  by  20  percent.  For  the  first  design  with  C/M  -  3.0, 
the  actual  C/M  ratio  to  be  used  becomes  3.6. 

From  example  1,  M  »  370.4  g  .  Therefore,  C  *  1533  g. 
The  height  of  the  Comp  B  disk  would  thus  be  6  cm. 

The  steel  tamper  weight  would  be  N  -  7M  *  2593  g 

The  height  of  the  steel  disk  would  therefore  be  2.54  cm. 

The  warhead  would  have  the  following  characteristics: 


o 


L/D  -  0.756 
C/M  -  3.6 
N/M  «  7.0 


(see  Fig  24b) 


If  C/M  *4.0  is  chosen,  the  warhead  characteristics 


are: 


L/D  -  0.779 
C/M  -  4.8 
N/M  -  2.5 

Height  of  Comp  B  disk  *  8  cm 

Height  of  steel  disk  -  0.943  cm  (see  Fig  24c) 

Of  the  two  designs,  the  first  one  would  be  preferable 
because  less  Comp  B  is  required,  given  that  the  cost  of  Comp 
B  is  far  greater  than  steel.  Compared  to  design  1,  it  is 
seen  that  the  addition  of  a  tamper  reduces  the  amount  of 
explosives  required  to  produce  the  same  velocities.  Also, 
as  the  tamper  gets  heavier,  the  amount  of  explosive  required 
gets  less. 

Since  no  experimental  data  are  available  for  this  type 
of  warhead,  it  is  impossible  to  determine  how  accurately 
the  Gurney  equation  estimates  the  velocity  distribution.  It 
is  assumed  that  the  method  provides  satisfactory  estimates 
for  preliminary  design. 

Example  3 

The  third  example  considered  is  similar  to  the  one  in 
Ref  2.  The  end  projector  is  treated  as  two  asymmetric 
sandwiches,  as  is  done  in  Fig  26.  It  is  assumed  that  the 
upper  part  will  be  given  a  velocity  of  approximately  100  m/s 
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by  the  Detasheet  layer.  It  is  necessary  to  find  the  charac¬ 
teristics  of  the  upper  part  which  will  give  a  velocity  of 
2900  m/s,  then  those  for  the  bottom  part  can  be  found.  The 
upper  part  is  made  up  of  a  layer  of  32  steel  cubes  surrounded 
by  a  plastic  ring,  a  layer  of  Comp  B  detonated  simultaneously 
around  the  perimeter  by  Detasheet,  1/16  inch  thick,  and  a 
steel  disk  beneath  the  Comp  B.  Since  this  is  similar  to 
example  2,  C/M  -  3.0  is  chosen.  This  implies  that  N/M  ■  7.0. 
Although  two  different  explosives  are  involved,  Con?)  B  and 
Detasheet,  the  amount  of  Detasheet  is  small  compared  to  the 
amount  of  Comp  B  and  the  Gurney  energy  for  Con?)  B  is  used  in 
the  computations. 

It  was  found  from  experimental  data  that  for  this  con¬ 
figuration  the  velocity  is  increased  due  to  the  amplifica¬ 
tion  of  the  detonation  wave  (Ref  19)  and  that  little  or  no 
correction  is  required,  even  though  pre- formed  fragments  are 
used.  Therefore,  in  this  case,  no  correction  is  used. 

For  the  upper  sandwich: 

Since  C/M  »  3.0  ,  C  =  1111.2  g  . 

Also,  N  -  2593  g  . 

The  steel  tamper  has  a  12.38  cm  diameter  and  is  2.78 
cm  thick. 

The  height  of  the  Detasheet  and  Comp  B  cylinder  is  5.28  cm. 
For  the  lower  sandwich  (Detasheet  between  steel  tan?)er  and 
upper  part) : 

The  weight  of  the  upper  sandwich  is  4074.6  g. 

The  weight  of  the  layer  of  Detasheet  between  the  upper 
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sandwich  and  the  lower  tamper  Is  approximately  33  g.  In 
this  case,  C/M  -  0.0081  .  To  obtain  a  velocity  of  100  m/s 

for  the  upper  part,  a  value  of  N/M  *  .3  is  chosen  because 
It  will  give  a  velocity  higher  than  100  m/s  and  interpolation 
in  the  tables  in  Ref  6  is  easier.  Also,  because  the  value 
of  C/M  is  so  small,  the  answers  obtained  are  highly  approxi¬ 
mate,  since  the  value  is  outside  the  bounds  of  the  Gurney 
theory.  However,  since  the  velocity  imparted  by  the  lower  part 
will  be  very  small  compared  to  that  from  the  upper  part,  the 
values  obtained  for  C/M,  N/M,  and  the  velocity  are  used  to 
get  an  estimate  of  the  size  of  the  bottom  tamper  plate.  By 
using  a  C/M  which  falls  within  the  range  of  applicability  of 
the  Gurney  equations,  more  accurate  answers  would  be  obtained, 
but  the  configuration  would  be  different  from  the  present  one 
if  a  fragment  velocity  of  3000  m/s  is  required.  This  is 
not  done  for  this  example.  Therefore,  the  warhead,  as  designed, 
would  have  a  diameter  of  12.7  cm  and  a  L/D  of  0.75  (see  Fig  24d) . 

Example  4 

Another  example  which  is  of  interest  is  one  where  there 
are  two  layers  of  fragments,  with  an  explosive  sandwiched 
between  the  fragment  layers .  This  warhead ,  producing  two 
fragment  clouds,  would  have  a  higher  probability  of  kill. 
Assuming  a  head-on,  or  beam  attack,  the  first  layer  would 
penetrate,  or  weaken,  the  target  while  the  second  layer 
would  destroy  it.  The  configuration  is  the  same  as  example 
3,  but  an  extra  layer  of  steel  fragments  enclosed  in  a 
plastic  ring,  backed  with  a  1/16  inch  layer  of  Detasheet, 
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are  added  to  the  warhead.  If  the  same  amounts  and  dimensions 
as  those  used  in  the  third  design  are  used,  the  fragment 
velocities  can  be  calculated  using  Gurney's  equation  for 
an  asymmetric  sandwich  (Eq  20) .  The  warhead  can  be  split 
into  three  different  sandwiches  to  calculate  the  velocities. 

The  first  part  would  be  made  up  of  the  steel  tamper  on 
the  bottom,  a  Datasheet  layer,  and  the  upper  part  comprising 
a  steel  disk.  Comp  B,  and  the  layers  of  fragments.  From 
example  3,  C  -  33  g  and  N  -  1222  g 

M  '  Ms  +  Mcomp  B  +  Mfl  +  Mf2  +  MDetasheet 


where 

Ms  -  steel  disk  -  2593  g 

Hcomp  B  "  mi- 2  g 

■  fragment  layer  *  359.8  g 

Mf2  »  uppermost  fragment  layer  *  370.4  g 

^Detasheet  ^  ^ 

Therefore,  M  «*  4476.4  g  and  consequently,  C/M  »  0.0074  , 

and  N/M  *  0.274  .  From  the  tables  in  Ref  6,  the  velocity  of 

the  upper  part  is  found  to  be  32  m/s.  Again,  the  C/M  falls 
outside  the  range  of  the  Gurney  equations,  and  the  velocity 
obtained  is  approximate.  Since  no  better  estimate  can  be 
obtained,  it  is  used  for  the  purpose  of  this  example. 

The  second  part  is  made  up  of  the  steel  disk,  Comp  B 
and  two  layers  of  fragments.  From  example  3,  C  =  1111.2  g  , 

N  -  2593  g,  and  M  ■  763.2  g  Therefore,  C/M  -  1.46  , 

^  and  N/M  “3.4  Again,  from  the  tables  in  Ref  6,  the  velocity 
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of  the  fragment  layers  is  2309  m/s .  When  the  initial 
velocity  imparted  by  the  bottom  layer  of  Detasheet  is  added 
to  the  fragments,  the  velocity  is  actually  2391  m/s. 

The  last  part  is  made  up  of  the  two  fragment  layers  and 
the  Detasheet  layer  between  them.  For  this  section,  C  *  33  g, 
H  *  370.4  g,  and  N  «=  359.8  g  .  Therefore,  C/M  *  0.089 
and  N/M  =  0.971  .  The  velocity  of  the  uppermost  fragments 

is  497  m/s.  For  the  bottom  layer,  it  is  514  m/s,  but  in  the 
opposite  direction  from  that  of  the  upper  fragments .  When 
the  velocities  that  have  already  been  imparted  to  the  frag¬ 
ments  are  added,  the  upper  layer  has  a  velocity  of  2888  m/s 
and  the  other  layer  has  a  velocity  of  1877  m/s.  These  are 
the  velocities  of  the  fastest  fragments;  those  at  the  edge 
of  the  projector  would  be  affected  by  the  end  effects  and 
would  have  lower  velocities. 

The  initial  velocity  of  3000  m/s  has  not  been  reached; 
the  best  way  to  attain  that  velocity  would  be  by  increasing 
the  amount  of  Coup  B  in  the  middle  layer  to  1435  g.  This 
would  give  the  top  fragment  layer  a  velocity  of  2500  m/s 
and  the  required  initial  velocity  would  be  reached  when  the 
other  velocities  are  added.  The  design  would  have  a  length 
to  diameter  ratio  of  1.03  (see  Fig  24e) . 

These  four  conceptual  designs  were  used  to  show  how 
the  modified  Gurney  method  can  be  used  in  preliminary  design. 
Experimental  testing  or  computer  simulation  would  be  necessary 
before  a  final  design  would  be  accepted. 
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VI  Conclusions  and  Recommendations 

It  is  possible  to  obtain  more  accurate  estimates  of 
fragment  velocity  distributions  by  using  modified  Gurney 
equations.  Basing  these  modifications  on  available  experi¬ 
mental  data,  it  is  concluded  that: 

1-  the  influence  of  end  effects  on  cylindrical  charges  with 
L/D<2  can  be  modelled  adequately  using  local  C/M's  which 
are  assumed  to  be  reduced  from  the  actual  C/M.  It  is 
assumed  that  a  cone  of  explosive  of  diameter,  D  (obtained 
from  the  value  of  L/D) ,  and  of  height,  2R,  is  removed  from 
the  initiating  end  and  that  a  cone  of  diameter,  D,  and 

of  height,  R,  is  removed  from  the  other  end, 

2-  the  influence  of  end  effects  on  end  projector  type  charges 
can  be  adequately  modelled  by  assuming  that  the  amount 

of  explosive  contained  is  a  cone  of  diameter,  D,  and  of 
height,  R  (obtained  from  the  value  of  L/D), 

3-  for  charges  initiated  at  one  end,  the  presence  of  pre¬ 
formed  fragments  can  be  modelled  by  computing  the  velocity 
distribution  using  a  C/M  which  is  80  percent  of  the 
actual  C/M, 

4-  for  end  projector  type  charges  where  the  explosive  is 
initiated  simultaneously  around  the  perimeter  of  the 
explosive,  no  reduction  in  the  C/M  is  required  to  correct 
for  the  presence  of  pre-formed  fragments, 

5-  and  that  use  of  these  modified  equations  will  yield 
reasonably  accurate  results  for  estimates  of  the  velocity 
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distribution.  The  modified  Taylor  equation  also  gives 
better  estimates  for  fragment  projection  angles  than  the 
usual  Taylor  equation. 

From  the  information  gained  in  this  project,  it  is 
recommended  that : 

1-  more  experimental  data  be  obtained  for  a  greater  variety 
of  L/D  and  C/M  for  cylindrical  and  end  projector  warheads. 
This  data  will  serve  to  improve  the  accuracy  of  the 
modifications  or  cause  them  to  be  changed, 

2-  the  modified  Gurney  equations  be  used  for  preliminary 
warhead  design, 

3-  and  that  final  design  be  based  on  experimental  testing 
or  computer  simulation. 
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Appendix  A 

Least  Square  Fit  for  Gurney  Energy ,  /7Z 

Each  explosive  has  a  particular  value  of  Gurney 
energy.  Also,  each  explosive  has  its  own  detonation  rate, 
dependent  on  density  of  packing  and  other  factors.  Since 
the  energy  appears  to  vary  linearly  with  the  detonation 
velocity,  a  least  square  fit  (Ref  6)  was  made. 

The  equation  for  the  straight  line  is: 

/2T  -  0.52  +  0.28  D 

where  D  is  the  detonation  velocity  of  the  explosive. 

The  equation,  along  with  experimental  data  from  Ref  4 
and  Ref  10,  are  plotted  in  Fig  25.  The  graph  abscissa  starts 
at  6.6  mm/ys  because  explosives  with  lower  detonation  rates 
are  of  little  military  use.  From  this  graph,  an  approximate 
value  of  the  Gurney  energy  can  be  found  for  any  explosive, 
given  its  detonation  velocity. 
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Appendix  B 


Sample  Calculation  of  Fragment  Velocity  for  End  Projector 


For  Shot  43  (Ref  2) 

C-4:  Weight 

Detonation  velocity 
Gurney  constant,  /7E 

Detasheet:  Weight 

Detonation  velocity 
Gurney  constant, 


"  284  g 

-  8.04  mm/us  (Ref  17) 

-2.79  mm/us,  from 
Fig  25 

-  54  g 

■  6.8  mm/us  (Ref  17) 

-  2.44  mm/us,  from 

Fig  25 


The  system  is  considered  as  two  asymmetric  sandwiches. 

The  first  is  comprised  of  the  plywood  disk,  Detasheet,  and 
all  that  is  above  (see  Fig  26a) ,  When  the  Detasheet  detonates , 
It  imparts  a  velocity  to  the  plywood  and  the  upper  sandwich. 
This  upper  sandwich  is  made  up  of  the  steel  disk,  C-4,  and  32 
Aluminum  cubes  imbeded  in  a  plastic  ring.  When  the  C-4  deto¬ 
nates,  it  imparts  a  velocity  to  the  steel  disk,  and  the 
plastic  and  fragments.  The  projector  has  a  12.7  cm  diameter. 

The  amount  of  Detasheet,  C^,  in  the  first  sandwich 
is  found: 


total  area  of  Detasheet  «  area  of  disk  +  area  around  sides 

-  126.71  +  76 

-  202.71  cm2 

CD1  -  54  x  202 *'7l  "  33*75  8  "  charge  wt  of  Detasheet 

The  amount  of  Detasheet  on  the  sides,  C^.  is  20.25  g. 


The  weight  of  the  steel  disk  Is  found: 


Density  of  steel,  pg  -  7.75  g/cc  (Ref  16) 

M,  -  p,  Vo!,  -  p,  -  624  g 

The  mass  of  the  A1  fragments  is  ^  -  32  x  2.26  -  72.32  g 

ttD  2 

Area  of  plastic  "  — -  area  of  cubes  **  95.74  cm  . 

Density  of  plastic,  pp  -  1.74  g/cc  (Ref  2) 

..Mp  -  Pp  Volp  -  158.79  g 

Therefore,  C/M  for  the  first  sandwich  is 

M  “  ”01.2  +  Ms  +  “c-4  +  *p  +  *M 
-  1159.4  g 
Cjj/M  -  0.0291 

Density  of  plywood,  ppl  -  513  g/cc  (Ref  16) 

Weight  of  plywood,  Mpl  -  ppl  Volpl  -  61.91  g 
Therefore,  Npl/M  "  "~lf  5$V4' '  “  0.0534 

From  the  tables  in  Ref  6,  the  speed  of  the  upper  part  is 
found  to  be 

VM  -  /2E  (0.041)  -  93  m/s 
For  the  upper  sandwich  (see  Fig  26b) 

C  *  CC-4  +  °D2  "  30425  S 

C/M  -  .488 
s 

Mass  of  plastic  and  fragments  ■  M  +  M^ 


c 


r 


-  231.11  g 

For  this  sandwich,  N  -  231.11  g  and, 

N/M  -  0.3704 

9 

Again  from  the  tables  In  Ref  6,  the  velocity  of  the 
fragments  Is 

VM  -  /7E  (.816)  -  2277  m/s 

But  since  the  Detasheet  gave  it  an  additional  velocity 
in  the  same  direction,  the  total  fragment  velocity  is 

VM  -  2277  +  98  «  2375  m/s 

Similar  calculations  can  be  made  for  Shot  No  66.  The 
computations  give  an  answer  for  of  2870  m/s. 
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Fig  26.  Cross  Section  of  Asynmetric  Sandwich 
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Conversion  Factors 


cm 

X 

.3937 

- 

inch 

g 

X 

.00022 

- 

pound  (av.' 

g/cc 

X 

.03613 

- 

lb/cu  inch 

mm/ vs 

X 

1000 

m 

m  /sec 

m/s 

X 

3.281 

m 

feet/sec 

Joseph  Yves  Charron  was  bom  on  26  August  1947  in 
Hull,  Province  of  Quebec,  Canada.  In  September  1966,  he 
was  accepted  to  the  Royal  Military  College  of  Canada, 

Kingston,  Ontario.  He  graduated  in  May  1970  with  the 
Degree  of  Mechanical  Engineering  and  was  granted  a  commission 
in  the  Canadian  Armed  Forces.  He  entered  the  Graduate 
Aeronautical  Study  Program  at  the  Air  Force  Institute  of 
Technology,  Wright-Patterson  AFB,  OH,  in  January  1978. 

Permanent  address :  DPCO/AERE 

National  Defense 
Headquarters 
Ottawa,  Ontario 


•T* 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  r*Ti«n  (>«•  EnlrriA) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

REPORT  number 

AFIT/GAE/AA/79S-1  * 

2.  GOVT  ACCESSION  NO- 

»  RECIPIENT'S  CATALOG  NUMBER 

a.  title  r«n<  Submit) 

Estimation  of  Velocity  Distribution 
of  Fragmenting  Warheads  Using  a 
Modified  Gurney  Method 

4.  TYPE  OF  REPORT  4  PERIOD  COVCREO 

*.  PERFORMING  oro.  REPORT  NUMBER 

T.  AUTHOR**) 

Yves  J.  Charron,  Capt,  CAF 

1.  CONTRACT  OR  GRANT  NUMBER*  > 

t  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Air  Force  Institute  of  Technology*" 
Wright-Patterson  AFB,  OK,  45433 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  0  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AFIT/EN,  Wright-Patterson  AFB, OH 

12.  REPORT  OATE 

11.  NUMBER  OF  PAGES 

14-  MONITORING  AGENCY  NAME  A  AOORESS*il  dllloranl  Iran  Control  lint  Off!  c) 

IS.  SECURITY  CLASS,  (ot  ihlm  topari) 

IS..  DECLASSIFICATION.  DOWNGRADING 
SCHEDULE 

It.  DISTRIBUTION  STATEMENT  foft/lti  Report) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ot  (he  abafracf  i%i »r*<f  In  Block  20,  U  dllforonl  from  Ropori) 

.  bS )  'x ... 

|  ptwctor  of  InformatiOMT  -43 

’  1  OCT  1 

1*.  SUPPLEMENTARY  NOTES 

19.  KEY  BORDI  fCanffrtui  an  rtvsns  aids  If  necessary  and  Identity  by  btack  number) 

End  Effects  Velocity  Distribution 

End  Projector  Warheads 

Fragments  Taylor  Formula 

-'•gurney 

SO.  * SSTRAC T^%^illnu.  on  r«v*ra«  ,tdm  If  nocommmry  mid  Identify  by  bio ek  nu mbit) 

The  Gurney  energy  method,  an  analytical  method  with 
which  to  compute  fragment  velocity  distributions  for  conti¬ 
nuous  casing  warheads,  is  reviewed  and  a  method  to  model 
velocity  losses  due  to  end  effects  on  cylindrical  charges 
and  end  projector  type  warheads  is  presented.  The  Taylor 
formula,  which  estimates  the  angle  of  projection  of  the 
fragments,  and  a  modified  Taylor  formula,  which  gives— ^ 

DD  1  jAN*?!  1473  *OtT|ON  or  1  NOV  81  IS  OBSOLETE  UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGE  Dmtm  Entoro*) 


SECURITY  CLASSIFICATION  QW  THIS  PACK(Wh*t  Dmim  Bntmd) 


*+■ 

— -^better  estimates,  are  also  presented.  The  method  is 

extended  for  use  with  casings  made  of  pre-formed  fragments. 
Finally,  four  conceptual  preliminary  designs,  based  on 
information  acquired  in  the  project,  are  investigated.^ 

v  The  modifications  made  to  the  equations  provide 
improved  results  and  the  examples  confirm  that  the 
modified  Gurney  method  is  a  quick,  inexpensive  tool  for 
use  in  preliminary  warhead  design.' 


S«CU*ITY  CLASSIFICATION  op  tu''  PACefNhFn  Dam  Cntanii) 


